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Abstract. We address the problem of energy dispersion of radiation pressure 
accelerated (RPA) ion beams emerging from a thin (solid) target. Two different 
acceleration schemes, namely phase-stable acceleration and multi-stage acceleration, 
are considered by means of analytical modelling and one-dimensional particle-in-cell 
simulations. Our investigations offer a deeper understanding of RPA and allow us to 
derive some guidelines for generating monoenergetic ion beams. 
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1. Introduction 

Interaction of ultra-intense laser pulses with thin foils offers interesting possibilities to 
generate energetic charged particles. So-called radiation pressure acceleration (RPA) of 
ion bunches has recently attracted a lot of interest as it may provide an efficient way to 
generate intense quasi-monoenergetic ion beams. In contrast to target normal sheath 
acceleration (TNSA) [1, 2], where ions are accelerated from the target rear surface (the 
front surface being the one irradiated by the laser pulse) in the electrostatic field built up 
by the laser-created hot electrons, RPA of ion beams relies on the efficient momentum 
transfer from laser photons to ions in a thin dense target, which reflects the incident 
laser pulse. RPA turns out to be a very efficient way to accelerate quasi-neutral ion- 
electron bunches up to potentially relativistic velocities while keeping energy dispersion 
small. 

The idea of accelerating (macroscopic) objects by use of laser radiation pressure 
was initially discussed by Marx [3] as a possible path toward interstellar space travel. 
Its application to efficient ion acceleration was first proposed in Ref. [4], where the 
authors show that, in order to observe efficient RPA with linearly polarized light, ultra- 
high laser intensities (II > 10 23 W/cm 2 ) are required. For lower intensities indeed, the 
strong laser-induced electron heating makes target expansion the dominant acceleration 
mechanism. In 2005, Macchi et al. demonstrated that using circularly polarized (CP) 
laser light strongly reduces electron heating thus allowing RPA to operate efficiently 
at lower intensities (II > 10 20 W/cm 2 ) [5]. After that, different teams have discussed 
independently the possibility to create quasi-monoenergetic ion beams by irradiating a 
thin target with an intense CP laser pulse [6, 7, 8, 9]. Many studies have followed, some 
proposing and/or revisiting different acceleration models or discussing the optimal laser- 
target parameters through one-dimensional (ID) and two-dimensional (2D) particle-in- 
cell (PIC) simulations [10, 11, 12, 13]. Multi-dimensional effects on the stability of the 
accelerated foil and their potential capability to improve RPA have also been discussed 
in Refs. [14] and [15]. Finally, first experimental confirmation of RPA has recently been 
claimed [16]. 

While the basic mechanisms of RPA are now well understood, the control of the 
energy dispersion in RPA ion beams has not been fully addressed. This paper is intended 
to provide the reader with guidelines how to achieve this control. To do so, we first 
recall the basic modelling of RPA of a thin foil and provide scaling laws concerning the 
maximum ion energy that can be reached as a function of the laser intensity or power 
(Sec. 2). Beyond this simple (macroscopic) modelling, we discuss the details of RPA 
of a thin foil as following from two complementary processes. For the thinnest foils, 
acceleration proceeds in the so-called phase-stable way [9]. For thick enough targets, 
acceleration occurs as a multi-stage process [6, 7]. In Sec. 3, we propose some refined 
models for both acceleration processes and extract the main requirements for controlling 
the resulting ion energy dispersion. These analytical findings are compared to numerical 
PIC simulations in Sec. 4. Finally, we present our conclusions in Sec. 5. 
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2. Basic Modelling 

In this paper, all quantities are normalized to laser and electron related units. Times and 
distances are normalized to the incident laser frequency ojl and wave-number = ul/c, 
respectively, and velocities are normalized to the light velocity c. Electric charges and 
masses are normalized to the electron charge e and mass m e , respectively. Densities 
are normalized to the critical density at the considered laser wavelength = l-nfki;. 
n c = eo m e ul/e 2 , where eo is the permittivity of free space. Electric fields are normalized 
to the Compton field Eq = m e cuL/e. Furthermore, we consider a CP laser pulse and 
introduce the incident laser pulse vector potential: 

A L (t, x) = cos(t — x) y + sin(t — x) z , (1) 

where y and z denote the two directions transverse to the laser propagation direction 
x. 



2.1. Macroscopic approach: the light-sail model 

A straightforward and elegant way to derive the energy ion gain during the acceleration 
of a thin target by laser radiation pressure is to consider the accelerated layer as a 
quasi- neutral light sail reflecting the incident laser pulse [1]. Acceleration then follows 
from momentum transfer from the laser photons to the ions. Assuming that all ions in 
the target have the same velocity vf in the laboratory-frame, the equation of motion 
for the target is obtained by equating the photon momentum flux to the ion momentum 
flux, which follows from the target acceleration (niodo) dpf^ /dt, where pf is the ion 
momentum in units of m e c, n i0 is the initial target ion density and do its thickness. For 
arbitrary target velocities, two effects - the reduction of the photon flux on the target 
and the Doppler shift lowering of the photon momenta in the target-frame - must be 
taken into account. Considering total reflection of the laser pulse in the target-frame 
one obtains: 

Ko do) j/P = 4[t - x t (t)] 7i 2 (l - v?) 2 , (2) 

where a\{t) is the normalized laser intensity, Xi(t) is the time-dependent position of the 
target moving with the velocity vf\ and 7$ = (1 — t^ 2 ) -1 / 2 is the associated Lorentz 
factor. The solution of Eq. (2) has been derived in Ref. [8] for a laser pulse with a step- 
like temporal profile [a\(t) = a 2 , for t > and a\(t) = otherwise], with the maximum 
laser field amplitude a : 

pf\t) = vrii [sinh0 — (4 sinh^)" 1 ] , (3) 

where = (1/3) sinh _1 (3 a\ tj (n^o mi do) + 2) and mi is the ion mass. 

At this point, we stress that, in the derivation of Eq. (2), (i) the electron momenta 
have been neglected, and (ii) the target is assumed to be thin enough to be accelerated as 
a whole, quasi-neutral bunch, but thick enough to support the laser radiation pressure. 
Both assumptions are discussed in more details in Sec. 3. 
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2.2. Scaling laws for the ion energy 

Let us now derive some scaling laws for the ion energy £, = (7; — 1) m, (in units of 
m e c 2 ). First, Eq. (3) can be simplified in the limit of non-relativistic ions (pf' 1), 
leading to pf\t) ~ alt/(n i0 d ) and ion energies: 

777; / Oq t 



This result suggests that the ion energy scales with the square of the laser intensity 
II — a o (i n units of c 3 m e n c /2), and more precisely with the square of the laser fluence 
(f>L = Jq a 2 L (t) dt [in units of c 3 m e n c / (2ujl)\. However, this scaling applies only for 
sufficiently short laser pulses. Indeed, if the laser pulse is long enough for the target 
to travel over a distance larger than the laser Rayleigh length Lr ~ w\ {wl is the 
transverse size of the laser focal spot), diffraction of the laser pulse must be accounted 
for. It sets in after a time t ~ y/2rrii do wl/clq, so that the final ion energy is limited to 
Si ~ <2q w 2 L /{nio do). For a sufficiently long laser pulse, the ion energy thus scales (only) 
linearly with the laser intensity. More precisely, one can introduce the normalized laser 
power onto the target, Pi ~ o>o w l [in units of m e n c c 3 /(2 fcfj], and we obtain that the 
final ion energy scales linearly with the laser power. 

These estimates suggest that relativistic ion energies can be reached for laser 
powers Pl > mi n iQ do. Considering characteristic diamond-like carbon (DLC) targets 
(mi = 12 x 1836, ni = 60 at Al = 1 fim) with thickness d/\i = 10~ 2 [17], this 
corresponds to a power P L > 3.5 x 10 4 , typically ~ 30 TW. Hence, current laser facilities 
could be considered for relativistic ion generation. 

Equation (3) can also be simplified in the limit of ultra- relativistic ions pf 1: 



2 \min i0 doJ 

The ion energy thus increases as t 1 ^ 3 . This characteristic evolution has been first 
observed in the original paper by Esirkepov et al. [4]. In the ultra- relativistic regime, 
~ 1, and diffraction sets in after a time t ~ w\. The final ion energy then scales as 
the power 1/3 of both the laser intensity and power. This is mainly because of photon 
red-shifting and the reduced photon flux onto the target due to its relativistic recoil. 



2 + \ 1/3 



3. Two RPA regimes 

A deeper insight into RPA of thin targets requires to investigate more closely the 
structure of the accelerated target. When exposed to an intense laser pulse, the target 
electrons are pushed forward into the target by the laser ponderomotive force thus 
forming a compressed electron layer (CEL) at the laser front. The formation of this 
layer occurs on a very short time as it involves only electron motion. Its characteristic 
position l c before ion motion sets in can be easily derived from equating the electrostatic 
pressure (Z n i0 l c ) 2 /2 to the radiation pressure on the CEL. Obviously, maintaining 
the target integrity requires l c to be smaller than the target thickness do. It is thus quite 
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natural to introduce the normalized parameter £ = Z d$ / (\/2 a$) . For £ < 1, the 
radiation pressure is so strong that it cannot be balanced by the electrostatic pressure 
inside the target. All electrons are expelled from the target which than undergoes 
Coulomb explosion. This regime of interaction has been studied in Refs. [19, 20, 21], 
where authors have considered its applications to both electron or ion acceleration. It 
is however not suitable for efficient RPA of thin foils, which requires £ > 1. 

In what follows, we analyze two regimes of RPA of thin foils for parameter values 
£ ~ 1 or £ > 1. 

3.1. Phase-stable acceleration of thin targets 

In the regime where £ ~ 1, all electrons are piled- up at the rear-side of the target and 
ion acceleration proceeds in the so-called phase-stable regime [9] (also referred to as self- 
organized double-layer acceleration [11] or coherent acceleration of ions by laser [18]). In 
this specific regime, ion acceleration can be described by considering that all electrons 
form a compressed layer at the target rear-side. For the sake of simplicity, we assume 
that the electron density in the CEL is constant and homogeneous: n ec ~ Z riiodo/l e , 
where l e < do is the CEL thickness. At this point, we stress that simple considerations 
on the balance between electrostatic and radiation pressures do not allow us to derive 
this thickness l e . For £ = 1, l e would indeed shrink to 0, which is prevented by the 
electron pressure that is not included in our model. 

As ions initially located outside the CEL are accelerated in an electrostatic field, 
which increases linearly in space, they are, a priori, not of interest for efficient generation 
of quasi-monoenergetic ion bunches. We should thus focus our attention on the ions 
located in the CEL which undergo acceleration in the monotonically decreasing field: 

E x (t,x) = E (t)^f^, (6) 

where d(t) is the position of the target rear-side at time t, and the maximum electrostatic 
field Eo(t) can be derived from the equilibrium condition of the CEL (in the frame 
moving with the target rear-side): 

l -{Zn l0 l c )E,{t) = al^ (l-^ } ) 2 , (7) 

where i> 4 - is the velocity (in the laboratory-frame) of the CEL, i.e. of ions accelerated 
in a phase-stable way, and % = (1 — vf^ 2 )~ 1 / 2 is the associated Lorentz factor. From 
this we obtain: 

E (t) = v / 2a o7 - (I""®) 8 - (8) 

Note that in this regime, where electrons are piled-up at the rear-side of the target, the 
accelerating field E (t) does not depend on the CEL thickness l e . A similar feature was 
discussed in Ref. [12]. 

The governing equation for the mean ion momentum can be easily derived by 
considering that, in the phase-stable regime, ions are accelerated, in the average, by an 
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electrostatic field E (t)/2. We obtain: 

which is nothing but Eq. (2) derived in the macroscopic model for £ ~ 1. 

As for the ion motion around the mean velocity vf\ it can be described as in 
Ref. [9]. Denoting Xi{t) the position of an arbitrary ion in the CEL with respect to the 
center of the CEL and considering that all ions move with a velocity close to the mean 
velocity vf \ we obtain: 

d 2 ZE (t) 
dt 2 Xl rrii l e 7 t ? 

Therefore, if both Eo(t) and 7; vary slowly on a time scale where Q 2 = 

Z Eo(t)/(mil e , yf), the ions in the CEL have an harmonic motion around the mean 
velocity. From this, one can infer the dispersion in ion velocities of the accelerated 
bunch Avi ~ l e Q in the frame moving with the CEL. Correspondingly, the relative 
energy dispersion for non-relativistic ions scales as: 



^Xi--—r4xi. (io) 



— 1/2 

This scaling with £ i ensures, with the small CEL thickness ( e « 1, the quasi- 
monoenergetic feature of the accelerated ion beam. In the ultra-relativistic regime, 
one would obtain ASi/Si oc £ i 5 ^ 2 . However, and as it will be shown in numerical 
simulations (Sec. 4.2), the reduction of radiation pressure in the target frame associated 
with its relativistic recoil does not allow for this acceleration regime to be maintained 
at ultra-relativistic velocities. Instead, acceleration will more and more evolve like the 
multi-stage process discussed in the next Sec. 3.2. 

Phase-stable acceleration of thin target thus opens an interesting path toward 
the creation of energetic quasi- mono chromatic ion beams. Nevertheless, there is one 
restriction which was not mentioned in the original paper by Yan et al. [0] that we want 
to address now. In this specific regime of laser-target interaction, the electron bunch is 
compressed at the target rear-side. A large energy can thus be stored in the electrostatic 
field. Once the laser pulse is turned off, this energy is transferred back to the electrons, 
which start to quiver around the target, thus inducing its adiabatic expansion and in 
turn widening the ion energy spectrum. 

To estimate the importance of this effect, we derive the energy stored in the 
electrostatic field (see also [12]). Neglecting the contribution of the CEL due to its small 
thickness and considering that the electrostatic field varies linearly, E x (x) = E (t) x/d(t) 
for < x < d(t), the energy stored in the electrostatic field reads: 

d ^E 2 (x) , a 2 4 / 



£es(t) ~ J ^ dX = f Ti [} - Vf) d(t) . (12) 

For non-relativistic ion velocities this quantity simplifies to S es (t) ~ n i0 d Si(t)/3, 
suggesting that the energy stored in the electrostatic field is of the same order as the 



Energy dispersion in radiation pressure accelerated ion beams 



7 



total ion kinetic energy. For such ion velocities, one should therefore expect a broadening 
of the ion energy distribution once the laser is turned off. For ultra-relativistic velocities 
however, £ es is found to remain much smaller than the total kinetic ion energy, and the 
effect of adiabatic expansion on energy dispersion is negligible. 

These theoretical predictions are compared to numerical simulations in Sec. 4.2. 

3.2. Multi-stage RPA of thicker targets 

The acceleration mechanism for £ ^> 1 can be described as a multi-stage process [6, 7], 
where the target undergoes successive hole-boring (HB) processes. For the sake of clarity, 
we first present the multi-stage process in the case of non-relativistic ion velocities. Then 
we generalize the procedure to the relativistic case. 

3.2.1. Non-relativistic ion velocities In a first stage, ion acceleration follows from the 
laser-driven HB of the immobile target. The laser acts on the target ions as a piston, 
moving deeper into the target with the velocity v p o and reflecting an increasing number 
of ions [22, 23]. The piston velocity can be derived easily by considering the balance 
of radiation and electrostatic pressure in the frame comoving with the piston. If the 
laser field amplitude ao and the target density are constant, v p o will not change 
in time. For non-relativistic ion velocities, we obtain v p o = a / a/2 m 8 n iQ . During this 
stage, the ion velocity in the laboratory-frame ranges between vf = (corresponding 
to ions having not been picked up by the laser piston) and vf\ = 2t> p0 (corresponding 
to ions that have been reflected once by the laser piston). Index 1 in the ion velocity 
denotes the first acceleration stage. This first stage lasts up to the time tq = do/vpo 
when the piston reaches the initial position of the back of the target. Ideally, at the 
end of this stage, all ions of the target have been accelerated to the velocity vf\ = 2 v p0 
in the laboratory-frame and the whole target has been accelerated as a quasi-neutral 
bunch. At this point, we restrict ourselves to thin targets with thickness d <C v p ot p 
(with t p the laser pulse duration), a necessary condition for the multi-stage process of 
RPA. Targets with larger thickness will only undergo HB. 

To describe the second acceleration stage, i.e. for times t > To, we consider that 
the whole target is moving with the velocity = vf\ = 2 v p o in the laboratory-frame. 
Then, in the frame moving with the target, ion acceleration proceeds in a way similar 
to laser-driven HB. During this stage, ion velocities range between and 2t> p0 m the 
target-frame, which transforms to = 2v p q and vf\ — vf + ^ v po = ^ v po i n the 
laboratory-frame. 

Then, if t p > 2 r , a third acceleration stage starts during which the ion velocity 
in the laboratory frame ranges between 4 v p o and 6 v p q. For sufficiently long laser pulse, 
this multi-stage process goes on so that at the j th stage, the ion velocity ranges between 
vfj_i = 2 (J — 1) v p o and vfj = 2jv P Q. Correspondingly, the ion energy ranges between 
2 (j — l) 2 miVp and 2j 2 mjt>p . Generation of quasi-monoenergetic ion bunches thus 
requires the acceleration process to occur over many steps, i.e. over a time t ^> r . 



Energy dispersion in radiation pressure accelerated ion beams 



8 



If this condition is satisfied, the ion energy and energy dispersion at a time t ~ j r 
(j > 1) read: 



The first equation corresponds to the ion energy evolution expressed by Eq. (4), which 
was obtained using the macroscopic model in Sec. 2.1 in the limit of non-relativistic ion 
velocities. The second relation gives an information about the energy dispersion of the 
ion bunch. 

3.2.2. Iterative procedure for relativistic ion velocities To extend the multi-stage model 
to higher ion velocities, relativistic effects such as the radiation pressure diminution on 
the target and dilation of the characteristic stage duration due to the target relativistic 
recoil have to be accounted for. 

The initial stage of ion acceleration is once more similar to the laser-driven HB 
of the immobile target. Accounting for relativistic effects, the piston velocity is v p = 
v p0 /(l + v p0 ). During this stage, ion velocities thus range between (not yet reflected 
ions) and vf\ = 2v p /(l + v 2 ) (reflected ions). Correspondingly, we obtain the minimum 
and maximum ion energies during the first stage: £^f n = and £™ ax = (7^ — l)mj, 
where 7^ = (1 — ) _1//2 . This stage ends when the piston reaches the initial position 
of the target rear-side, i.e. after a time tP = do/v p . In contrast to the classical limit, 
however, ions do not reach exactly twice the piston velocity. As a consequence, the 
target thickness, and consequently the ion density, have changed. At the end of this 
first acceleration stage, we have d\ = vf^^ — d and n^i = n i0 do/d±. 

Let us now consider the j th stage of the acceleration process (j > 1). The 
target velocity (in the laboratory frame) is the velocity that ions have reached at 
the former j — 1 stage : v!jp = vfj_ 1 . In the frame moving with the target, 
the laser radiation pressure is therefore reduced by the factor 7 t 2 (1 — vf^) 2 [where 
7t = (1 — v® 2 ) -1 / 2 ]. The piston velocity in the target frame thus reads vjp = v'/(l +v'), 
where v' p = v p q (nioA^j-i) -1 ^ 2 It (1 — v 1 ) . It follows that the reflected ions have a 
velocity vf = 2v p ^ /(l + v p ^ 2 ) in the frame moving with the target, which transforms 
in vfj = (vf + v\ ; )/(l + v\ v\ ) in the laboratory- frame. The duration of the stage 
in the target frame is easily computed as rf 1 = dj-i/v p \ while one has to account 
for time dilation in the laboratory frame = "ft^- As for the target thickness and 
density at the end of the stage, they have to be recalculated in the target frame as 
dj = i> t - Ts — dj-i and n,j = n i0 d /dj. 

Following this procedure, we can compute the temporal evolution of the minimum 
and maximum ion energies for arbitrary values of the parameters v pQ and do. The 
comparison of this multi-stage model with the macroscopic model of Sec. 2.1 is given in 
Fig. 1 for several values of v p q. Predictions from the multi-stage model match perfectly 
with results from Sec. 2.1 for v p q = 0.01, where the characteristic ion energy evolution 
oc (t/r ) 2 is recovered (Fig. la). For higher values of v p0 (Fig. lb-d), a still good 



£i(t) = 2miDp (t/r ) 2 , 
A£i(t) = (2T Q /t)Si(t). 



(13) 
(14) 
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agreement is found between the two models. A small discrepancy can however be 
observed in the ion mean energy J, but it remains small compared to the predicted energy 
dispersion. It must also be noted that the characteristic duration of an acceleration 
stage, for relativistic ion velocities, can be strongly dilated in the laboratory-frame, 
which has an important effect on energy dispersion. These analytical predictions are 
compared to numerical simulations in Sec. 4.3. 

4. Numerical simulations 

Numerical simulations of RPA of thin targets have been performed using the PIC code 
PICLS [25]. In order to make a direct comparison with our analytical model, only 
one-dimensional in space and three-dimensional in velocities (1D3V) simulations are 
presented. This choice is also justified as it has been shown that a quasi-lD geometry 
is required to avoid strong electron heating and partial transparency of the foil, which 
may prevent quasi- monoenergetic ion beam generation [7]. Our study thus provides us 
with necessary, albeit not sufficient, conditions for creating monoenergetic ion beams. 

In our simulations, a circularly polarized laser pulse is focused at normal incidence 
on a thin, fully ionized, Carbon target with density = 25 and Z — 6. The target 
is located at a distance 2 from the left boundary of the simulation box (the laser 
propagates from left to right). Both the incident laser field amplitude and the target 
thickness are varied to explore different regimes of RPA in thin foils. 

4-1. Optimal target thickness 

Figure 2 shows the ion mean energy at an instant t ~ 10 tl after the beginning of 
the interaction for £ > 1. In the case £ < 1, the simulations indeed confirm that all 
electrons are removed from the target due to the strong radiation pressure, leading to 
the Coulomb explosion of the non-neutralized ion layer. For £ > 1, one observes that 
the ion mean energy increases as the target thickness decreases: the lighter the target, 
the higher ion energy one can reach. Such a result was already discussed in Ref. [10, 12] 
with the conclusion that the regime, £ ~ 1, is the optimum case for high energy ion 
generation. 

4-2. Ion acceleration in the phase-stable regime 

Let us have a more detailed look into ion acceleration in the phase-stable regime £ ~ 1. 
Figure 3 presents the time-resolved energy spectra obtained in simulations for different 
incident laser field amplitudes a = 10 — 100 and £ ~ 1 (the target thickness is adjusted 
for each laser amplitude). Theoretical predictions for the ion mean energy from the 
macroscopic model (Sec. 2.1) are superimposed on the numerical results. A rather 

| Our multi-stage model predicts a slightly higher ion energy than the macroscopic model. This is 
because, in our multi-stage model, we have a step-like decrease of the radiation pressure onto the 
target, while radiation pressure is continuously decreasing in the macroscopic model. 
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good agreement is found between theory and simulation. However, while for ao = 100 
the ion mean energy evolves exactly as predicted by the macroscopic model, numerical 
simulations at lower field amplitudes (ao = 10 — 20) show higher ion energies than 
estimated analytically. The reason of this discrepancy can be found in Fig. 4 and the 
corresponding movies, where details of the temporal evolution of the target structure 
during the acceleration process are presented for a = 10 and a = 100. For the lower 
laser field amplitude, a$ = 10, the CEL is not totally opaque to the laser field, which is 
partly transmitted (the foil transmittance is here T ~ 20 %). The ponderomotive force 
on the target rear-side is thus non zero and some electrons can escape from the target 
into the vacuum behind. This widens the electrostatic field distribution and increases its 
average value at the target rear-side thus leading to the observed increase of ion energy. 
In contrast, for ao = 100, the target is partly transparent only during a short time at 
the beginning of the interaction. Once the target has reached a relativistic velocity, 
the radiation pressure in the target frame is lower and electrons remain confined in the 
target. On such long times, the CEL does not actually stay at the target rear-side and 
ion acceleration becomes more similar to multi-stage acceleration. 

Figure 3 also reveals the quasi-monoenergetic ion distribution during the 
acceleration process. This is underlined in the different panels of Fig. 3 presenting 
the ion energy at the end of the interaction process (panels 3e-h), as well as in Fig. 5 
which shows the relative ion energy dispersion as a function of the ion mean energy for 
ao = 10 — 100. Figure 5 also confirms the characteristic dependency of the relative energy 

—1/2 

dispersion, which is proportional to S i as predicted by Eq. (11). Equation (11) also 
predicts that A£j/£j should scale as \/l e a , and one could naively expect for a given 
ion mean energy, that ASi/Si scales as the square- root of the laser field amplitude ao- 
Figure 5 however shows that the dependence on ao is stronger. This is because the CEL 
thickness l e actually increases with ao for a fixed value of £ oc do/ao- 

Let us now discuss the fraction fi of ions in the monoenergetic peak. A naive 
estimate can be derived from the semi-microscopic model presented in Sec. 3.1 by 
considering that only ions initially located in the CEL participate in the phase-stable 
acceleration so that fi ~ 1 — This estimate would suggest that only a very small 
fraction of the target ions participate in the acceleration process. Simulations however 
show quite the contrary: /j ranges between 0.45 and 0.67 for ao = 10 — 100. There are 
two reasons for such a high number of accelerated ions, (i) While in the model presented 
in Sec. 3.1 the CEL thickness shrinks to as £ — > 1, the electron pressure in the CEL 
actually increases during the compression by the laser pulse and prevents its collapse. 
As a result, the CEL thickness is much larger than expected in the model and so is the 
fraction of accelerated ions, (ii) Furthermore, some ions initially located outside of the 
CEL can still be injected (after some time) in the CEL and thereafter participate in the 
phase-stable acceleration. While creation of the CEL is almost instantaneous (because 
of electrons relativistic velocities, it occurs on a characteristic time ~ d /c <C 1), ions 
react to the strong electrostatic field on a longer time scale of the order of the inverse 
ion plasma frequency. A self-maintained structure made of both the CEL and ions 
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is formed and it is this structure which is accelerated in a phase-stable way. The 
correct modelling of this accelerating structure is very challenging as it would require 
to describe self-consistently the ion and electron motions. While this has been done for 
ion acceleration in the HB regime in Ref. [23] by developing a quasi-stationary model, 
this is particularly difficult under current conditions as the quasi-stationary hypothesis 
does not hold. What is also clearly highlighted in the two movies is that some ions, 
which are not initially located inside the accelerating structure, can be injected into it 
after some time. Ions located outside the CEL indeed 'see' a constant accelerating field 
and may actually catch up with the CEL after some time. For non relativistic ions, one 
can easily estimate that only ions initially located at a position Xj > d /2 can reach the 
CEL and their fraction /, cannot exceed 50 %. Accounting for relativistic effects allows 
for a larger fraction of reinjected ions (which increases in time, according to simulations, 
up to 67 % for a = 100). 

Finally, we want to point out that, once the laser is turned-off, the energy stored in 
the electrostatic field goes back to electrons which start quivering around the accelerated 
ions thus driving its adiabatic expansion. If the electrostatic energy is of the same order 
of magnitude as the total ion energy, this can strongly enhance the final ion energy 
dispersion. To investigate this effect in-depth, we have plotted the total energy of 
ions in the monoenergetic peak and the energy stored in the electrostatic field as a 
function of time for oq = 10 — 100 in Fig. 6. For small laser amplitudes ao = 10 — 20 
and correspondingly non-relativistic ion energies, a non-negligible fraction of the energy 
(~ 30 %) is stored in the electrostatic field. As a consequence, the corresponding ion 
energy spectra (Figs. 3i,j) are considerably wider at the end of the simulation. On the 
contrary, and as predicted in Sec. 3.1, this effect is negligible for sufficiently large ion 
energies. For ao = 40 — 100, most of the energy is stored in the accelerated ion bunch 
and no-enhanced ion energy dispersion is observed (Figs. 3k,l). 

In addition, we notice that this effect can also be mitigated at low ion energy by 
using a more sophisticated laser pulse temporal profile, e.g. by considering Gaussian or 
hyper-Gaussian pulses. Figure 7 shows the ion energy spectra obtained using a 6 t/l -order 
hyper-Gaussian or Gaussian laser pulse with similar maximum field amplitude ao = 10, 
fluence and full-width at half-maximum. The enhanced energy dispersion due to the 
long-time behavior of electrons is strongly mitigated as the pulse is slowly turned off. 

4-3. Ion acceleration in the multi-stage regime 

PIC simulations in the regime of multi-stage acceleration (£ 3> 1) are now discussed. 
Figure 8 shows the time resolved ion energy spectrum extracted from simulations with 
£ = 4 and different laser field amplitudes (ao = 10 — 100). In these simulations, the 
laser temporal profile is constant over a duration corresponding to 20 to ( we recall 
that to — do/ v po is the characteristic duration of an acceleration stage). Theoretical 
predictions from the model developed in Sec. 3.2 are superimposed to the numerical 
results. A rather good agreement on the minimum and maximum ion energies as a 
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function of time is found between theory and simulations: our model allows to correctly 
predict both the ion mean energy and energy dispersion. For cases with a$ = 40 and 
ao = 100, ions quickly reach relativistic energies and the dilation in the laboratory 
frame of the characteristic stage duration becomes obvious. Also, the duration of the 
laser-target interaction (and therefore of the acceleration process) increases as the foil 
reaches larger velocities (Fig. 8). 

Additional details on the target structure as well as the ion and electron phase- 
space distributions during the acceleration process are given in Figs. 9 and 10 and in the 
corresponding movies. The transition between successive acceleration stages is clearly 
visible in the movies. For the case ao = 10 for instance, the first HB stage terminates at 
t ~ 16 tl- As ions are reflected again and again by the laser piston, their distribution in 
phase-space becomes more and more complex in contrast to what was observed during 
phase-stable acceleration (compare for instance Figs. 9c, d and Figs. 4c, d). 

Furthermore, the ion energy spectra at the end of the simulations (Figs. 8i-l) are 
quite similar to those obtained at the end of the laser-target interaction (Figs. 8e-h). In 
contrast to phase-stable acceleration, the ion energy distribution here is not sensitive to 
the late time behavior of the electrons. This is due to the small fraction of laser energy 
which is stored in the electrostatic field at £ ^> 1. 

This excellent control of the ion beam spectrum as well as the large fraction of 
accelerated ions (/, exceeds 90 % in these simulations) make this regime of acceleration 
especially attractive for high-quality ion beam generation. Numerical simulations also 
suggest that RPA proceeds in the multi-stage regime as soon as £ > 2, which makes 
this robust mechanism more likely to be observed in experiments than phase-stable 
acceleration. 

Finally, as RPA in the multi-stage regime follows from successive HB of the 
target, one may suggest an additional source of energy dispersion for large laser field 
amplitudes and/or rather thick targets. Recent studies have indeed underlined the 
non-stationary behavior of the laser piston at high laser intensities giving rise to the so- 
called piston oscillations [23, 24]. This phenomenon implies large-amplitude (typically 
~ 30 %) oscillations of the maximum electrostatic field in the laser piston leading to an 
enhanced energy dispersion of ions reflected during the HB process. While the origin 
of these oscillations is yet not fully understood, some of their characteristic features are 
known: (i) their characteristic period is of the order of the inverse ion plasma frequency 
(ujpi = \J Z 2 nio I vrii in our normalized units), (ii) they appear after a characteristic 
time t s which is shorter for larger laser field amplitudes and/or target densities. For 
instance, PIC simulations of HB of a thick carbon foil with ion density = 25 by a 
CP laser pulse with ao = 40 indicate that oscillations in the electrostatic field with a 
characteristic period ~ 1.4 tl and ~ 35 % amplitude with respect to the maximum field 
strength occur after a time t s ~ 5 tl- Increasing the laser field amplitude to ao = 100 
does not change the relative amplitude of the oscillations or their period but shortens 
the characteristic time of their appearance to t s ~ 3.7 tl. 

One could therefore fear that these piston oscillations widen the energy spectrum 
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during the multi-stage acceleration. Clearly, this is not the case in the simulations 
presented in Fig. 8, where the energy dispersion is well described by our multi- 
stage model (which does not account for the piston oscillations). However, for these 
simulations, the characteristic duration of an acceleration stage to ~ ^o/^po — 6 tl is 
rather short and the piston oscillations have scarcely the time to develop. Therefore, 
we have performed a simulation with a thicker target (£ = 10 and ao = 40). Results 
from this simulation are presented in Fig. 11. During the first acceleration stage, for 
t = — 16 Tl, oscillations in the maximum electrostatic field are clearly visible after 
a time t ~ 5 T£ (Fig. lib) and result in a rather complex ion phase-space distribution 
(Fig. 11c as compared to Fig. 10a). Nevertheless, these regular oscillations disappear at 
the end of the first stage (t > 16 tl), which we attribute to the already perturbed target 
configuration in the second acceleration stage. There are still some large variations 
in the maximum value of the electrostatic field, but these variations are generic to all 
realistic simulations we have performed in the multi-stage regime. They follow more 
from the global ion dynamics in the target than from the piston oscillations themselves. 
Hence, these results suggest that the so-called piston oscillations are not a concern for 
the control of ion energy distributions during RPA in thin foils. 

5. Conclusion 

A detailed study of ion energy dispersion in RPA ion beams has been presented using 
both analytical modelling and 1D3V PIC simulations. The description proposed here 
allows for a greater insight in the details of RPA of thin foils than that available from 
the standard macroscopic light-sail model. In particular, it provides us with necessary 
conditions for quasi-monoenergetic ion beam generation. 

Two RPA regimes are identified depending on the dimensionless parameter £ which 
determines, for a given laser field amplitude, the target thickness. For both regimes, the 
models we have developed allow to recover the ion energy temporal evolution obtained 
by considering efficient momentum transfer from the laser photons to the target ions (the 
usual light-sail model). By accounting for the target structure during the acceleration 
process, we have gained a deeper insight into the ion energy dispersion. For £ ~ 1 (thin 
targets), RPA proceeds in the phase-stable regime introduced in Ref. [9]. Two sources 
of energy dispersion have been identified in this regime: the electric field inhomogeneity 
in the accelerating structure, and the adiabatic foil expansion due to the late time 
electron behavior. This later process is mainly important for low energy ions and its 
effect can be mitigated by using smooth temporal laser profiles. For £ > 2 (thicker 
targets), ion acceleration proceeds in the multi-stage regime originally discussed in 
Refs. [6] and [7], and for which we have developed a relativistic model. Ion energy 
dispersion in this regime is mainly determined by the number of acceleration stages. 
Small energy dispersion can thus be achieved by using long enough laser pulses. 

Hence, this work suggests that using moderately intense (and long) laser pulses is 
preferable for monoenergetic ion beam generation. The lower limit on the laser intensity 
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actually follows from the need to accelerate the ions to the desired energy before the 
target escapes from the laser focal volume. As for the optimal target thickness, the 
thinner the target is, the higher ion energies one can reach. However, maintaining 
the target integrity in the regime £ — > 1 might be experimentally difficult. Any 
nonuniformity in the laser intensity profile may lead to the removal of electrons and 
the resulting Coulomb explosion of the target. Also, Rayleigh- Taylor like instabilities, 
which have been observed in 2D simulations [7], may be more detrimental for very thin 
targets. Since simulations suggest that RPA occurs in the multi-stage regime as soon 
as £ > 2, phase-stable acceleration may thus be difficult to achieve in experiments. 
We therefore expect multi-stage acceleration to be the practically relevant acceleration 
mechanism. 
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Figure 1. Ion maximum and minimum energy predicted by the multi-stage model 
(dashed curves) and comparison to predictions from the macroscopic model (gray solid 
curves) for: a) v p o = 0.01, b) v p o = 0.1, c) v p q = 1 and d) v p q = 10. 
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Figure 2. Dependence of the ion mean energy on the normalized parameter £ 
(proportional to the target thickness) ~ 10 tl after the beginning of the interaction. 
The carbon target has density Zn^ = 150 and the laser field amplitude is clq = 10 
(o), ao = 20 (□) and ao = 40 (a). Red dashed lines correspond to predictions from 
the macroscopic (light-sail) model [solutions of Eq. (2)]. 
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Figure 3. a-d) Timc-rcsolvcd ion energy spectra, e-h) Snapshot of the energy 
spectrum at the end of the laser-target interaction, i-1) Snapshot of the energy 
spectrum at the end of the simulation. The laser field amplitude is: a,e,i) ao = 10, 
b,f,j) ao = 20, c,g,k) ao = 40 and d,h,l) ao = 100. Red dashed lines show theoretical 
predictions from the macroscopic model. 
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Figure 4. Details of the target 10 tl after the beginning of the interaction for 
a,c,e) a = 10, and b,d,f) ao — 100. a,b) The laser field amplitude is shown in 
magenta, the electrostatic field in green and the ion and electron densities in red and 
black, respectively. c,d) Ion distribution in phase-space. e,f) Electron distribution in 
phase-space. See also Movies 1 and 2 available from www.pks.mpg.de/~grech/RPADE/ 
mviel.avi [0.8 MB] and www.pks.mpg.de/~grech/RPADE/mvie2.avi [1.6 MB]. 
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Figure 5. Relative energy dispersion during the phase-stable acceleration as a function 
of the ion mean energy. For ao = 10 (o), a = 20 (□), a = 40 (A) and ao = 100 (V). 
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Figure 6. Temporal evolution of the total kinetic energy of ions in the monoenergctic 
peak (solid curves) and the energy stored in the electrostatic field (dashed curves) for: 
a) cto = 10, b) ao = 20, c) ao = 40 and d) ao = 100. Gray lines show theoretical 
predictions from Sec. 3.1. 




Figure 7. Same as Fig. 3 for £ = 1, ao = 10 and two different laser intensity profiles: 
a,b,c) for a 6 t?i -order hyper-Gaussian profile, d,e,f) for a Gaussian profile. Red dashed 
lines in panels a) and d) show theoretical predictions from the light-sail model [solutions 
of Eq. (2)]. Gray dashed lines in panels b,c,e,f) show the energy spectra obtained using 
a rectangular laser pulse profile. 
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Figure 8. Same as Fig. 3 for £ = 4. Red dashed lines in panels a-d) show theoretical 
predictions from the iterative model (Sec. 3.2). 
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Figure 9. Details of the target 20 tl after the beginning of the interaction for 
a,c,e) a = 10, and b,d,f) ao — 100. a,b) The laser field amplitude is shown in 
magenta, the electrostatic field in green and the ion and electron densities in red and 
black, respectively. c,d) Ion distribution in phase-space. e,f) Electron distribution in 
phase-space. See also Movies 3 and 4 available from www.pks.mpg.de/~grech/RPADE/ 
mvie3.avi [1.6 MB] and www.pks.mpg.de/~grech/RPADE/mvie4.avi [3.1 MB]. 
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Figure 10. Ion phase-space at different times after the beginning of interaction: a) 
3 tl, b) 10.6 tl, c) 18.2 T£, d) 25. 8tl and e) 33.4 T£. Each time corresponds to a 
different acceleration stage. In this simulation, ao = 10 and £ = 4.0. The color dots 
follow test ions during the acceleration process. The vertical dashed line shows the 
position of the laser piston (position of the maximum electrostatic field). 
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Figure 11. a) Time-resolved ion energy spectra. The dashed lines show theoretical 
prediction from the multi-stage model, b) Temporal evolution of the maximum 
electrostatic field. Vertical gray lines indicate the end of the different acceleration 
stages as predicted from the multistage model, c) Ion phase-space 14.6 tl after the 
beginning of the interaction. The vertical red line shows the position of the laser piston 
(where the electrostatic held is maximum). 



